Strict Physicality of Global Weak Solutions of a 
Navier-Stokes Q-tensor System with Singular 

Potential 



By Mark Wilkinson 



(N 

O 
(N 

> 
O 

^: 

(N 



CM 

< 



> 

oo 

o 



(N 



Abstract 

We study the existence, regularity and so-called 'strict physicality' of weak solu- 
tions of a coupled Navier-Stokes Q-tensor system which is proposed as a model for the 
incompressible flow of nematic liquid crystal materials. An important contribution to 
the dynamics comes from a singular potential introduced by Ball and Majumdar 
[3] which replaces the commonly employed Landau-de Gennes bulk potential. This is 
built into our model to ensure that a natural physical constraint on the eigenvalues 
of the Q-tensor order parameter is respected by the dynamics of the system. 

In this article, we construct global-in-time weak solutions to the following coupled 
Navier-Stokes Q-tensor system on the d-dimensional torus, which is an adaptation of a 
model of Beris and Edwards [4]: 

^(x, t) + (u(x, t) ■ V) Q{x, t) - S{Q{x, t),Vu{x, t)) 



r{LAQix,t)-e^iQix,t)) + ^tr 



dQ 



{Qix,t)) 



I + KQ{x,t)) , 



(S) 



du 
'dt 



{x, t) + {u{x, t) ■ V) u{x, t) = iyAu{x, t) — Vp{x, t) + div (r(x, t) + a{x, t)) 



V • u{x,t) = 0. 

The dimension d is either 2 or 3, Qisadxd matrix- valued map and T, L,9, KjL' > 
are constants. Moreover, ^p is a convex map defined on symmetric and traceless d x d 
matrices, whose construction was recently given by Ball and Majumdar [3] and which 
ensures the eigenvalues of the tensor field remain in a so-called physical interval The 
'tumbling' term S and the stress tensors r and a are defined in Section 2. Employing 
the map tp in system (S) allows us to infer that weak solutions Q{-,t) belong to L°° for 
t > 0. By virtue of the convexity ip, a suitable comparison principle argument allows 
us to infer strict physicality of weak solutions (see Section 4, Stage III below), which in 
turn allows us to prove higher regularity of global weak solutions in dimension 2. 



(^) Correspondence to: wilkinson@maths.ox.ac.uk. Address: Oxford Centre for Nonlinear Partial Differential 
Equations^ Mathematical Institute, Oxford University, 24-29 St Giles', Oxford, OXl 3LB. 
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We now discuss the modelling problem in condensed matter physics which motives 
our study of the above coupled system. The analysis of system (S) begins in Section 2, 
so the mathematically-minded reader may wish to skip Section 1 below. 

1. Motivation: Order Parameters and Physical Constraints 

Nematic liquid crystals form a class of condensed matter systems whose constituent 
rod-like molecules give rise to rich nonlinear phenomena, such as isotropic-nematic 
phase transitions. Thermotropic nematic liquid crystals form a subclass whose optical 
properties change dramatically with variation of system temperature. Above a certain 
temperature threshold the constituent molecules are randomly oriented (the isotropic 
phase), whereas below this threshold they tend to lie in locally preferred directions (the 
nematic phase). 

It is a formidable problem to formulate a mathematically rigorous model of such 
phenomena based on classical molecular dynamics in the continuum. This is, naturally, 
due to the high dimensionality of the phase space in which the dynamics take place. 
Due to the forbidding complexity of such models, one is led to consider more tractable 
mesoscopic or macroscopic models which are built by employing the general principles 
of thermodynamics and classical mechanics. 

1.1. Order Parameters. One of the first steps to take when formulating such 
effective static or dynamic continuum theories is to select an order parameter that 
captures the essential small scale structure in nematic systems. For meso-scale models 
such as the Maier-Saupe theory of statics [5] or dynamic Doi theories [9] , one typically 
posits as an order parameter a probability density p on the unit sphere S*^"^, which is 
to be regarded as encoding the average molecular orientation at material points in the 
spatial domain. As regards macro-scale models, there are a number of competing order 
parameters in the literature. Within the well-established Oseen-Frank [12] or Ericksen- 
Leslie theories [15], for example, one studies the respective static and dynamic properties 
of vector fields n : M*^ — > M'^ taking values in the unit sphere, where n(x) € S^~^ is known 
as a director. 

The director-field formalism may however be viewed as restrictive, as the order 
parameter n cannot account for hiaxiality of liquid crystal configurations. The Nobel 
prize-winning physicist Pierre-Gilles de Gennes suggested a traceless and symmetric 
dxd matrix with real components to be more appropriate in this regard for the modelling 
of nematic materials. 

1.2. The Q-Tensor Order Parameter. We now present some basic properties 
of the Q-tensor order parameter. For further details, consult De Gennes and Frost 
[7], Majumdar [18] or Newton and Mottram [19]. In what follows, the spatial di- 
mension d will be either 2 or 3, the two-dimensional case corresponding to, for instance, 
thin films of nematic materials. 
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1.3. Physical Constraints: an Eigenvalue Inequality. Suppose that to each 
point X in a material domain 17 C M'^ we associate a probability density function fx 
on molecular orientations which lie in S"^"^. In order to capture the Z2 'head-to-tail' 
symmetry of nematic molecules, each density is endowed with the antipodal symmetry 
fx{^) = fx{~^) for all uj G The Q-tensor order parameter Q{x) at x € is 

defined to be 

Qix) := [ (uj ^UJ - -l') fa:(uj)duj. 

It is the normalised matrix of second-order moments of the probability measure on S'^"^ 
with density fx- One may quickly check that Q{x) is a member of 

Symo(d) := {Q G M.'^'"^ : = Q and tr [Q] = o} . (1) 

The term —1/dI (which contains no information about the system) is included by 
convention in the definition so as to render the Q-tensor identically zero when fx is 
the uniform distribution on the unit sphere, corresponding to the isotropic phase of a 
nematic material. In this way, the Q-tensor order parameter should be interpreted as 
a crude measure of the deviation of a nematic system from isotropy. 

If one interprets a member of Symg(d) as a normalised matrix of second moments 
of some probability measure on nematic orientations, the eigenvalues of this order pa- 
rameter are then constrained in the following manner. By the spectral decomposition 
theorem, every matrix Q € SymQ(d) has the representation 

Q = Aiei (g) ei + ... + XdCd e^, 

where each is a unit-norm eigenvector of Q with corresponding eigenvalue Aj. If Q 
arises from some probability density /, we also have the identity 

J (^uj 0UI — i/^ f{uj) duj = XjCj (g) ej. 

Applying the action of the above matrices to a fixed eigenvector and then producing 
the result again with e^, one discovers the equality 

Xk= [ (w • ekff{uj) dhj - 3, 

from which one quickly deduces that is constrained to satisfy the inequality — l/ii < 
Afc < 1 — 1/d for k = 1, ...,d. The cases of equality correspond to perfect crystalline 
nematic alignment, and so are excluded. Thus, the eigenvalues Xi,...,Xd of any such 
Q £ Symg((i) should satisfy 

-^<Afc<l-i for k = l,...,d. (2) 

We subseqently refer to all Q-tensors whose eigenvalues satisfy inequality (2) as physical. 

From a modelling perspective, the physicality requirement (2) on the order param- 
eter Q presents a problem. It is often not clear (and sometimes not even the case) that 
solutions of static and dynamic theories of nematics which employ the Q-tensor order 
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parameter respect this physicality condition. Clearly, an issue then arises as to how one 
should interpret the solutions of such theories in a meaningful way. 

One well-studied theory that employs the Q-tensor order parameter is the static 
Landau-de Gennes theory. The Landau-de Gennes energy in the case d = 3 is given by 

E^,g[Q] ■■= (^IVQP + ^tr _ + £ (tr [Q^])') dx, 

for a suitably smooth domain 17 C and with a € M and b,c> 0. In [18], Majumdar 
shows that bulk energy minimisers in the spatially-homogeneous Landau-de Gennes 
theory (VQ = 0) do not always respect the above physical constraint on eigenvalues. In 
particular, when laboratory values for the liquid crystal material MBBA are strapped 
to the material-dependent constants (a, 6, c), one finds that (2) is not respected 2°C 
below the isotropic-nematic phase transition temperature. 

Aside from such issues related to the practical value of the unconstrained Landau- 
de Gennes theory for modelling real nematic materials, one can appreciate from this 
simple example of Majumdar that it is not immediately clear in the more complicated 
case of inhomogeneous spatial profiles whether or not minimisers of i?LdG subject to 
physical boundary conditions are physical pointwise on the spatial domain 17. 

Naively, one might hope to circumvent this problem in general simply by restricting 
the class of candidate minimising maps to 

A := |(5 G H^{^}) • ~^ < < 1 — - for A; = 1, d and almost every x € $7 

for instance. Even if one were to consider this new model, the problem remains for the 
case of dynamics. It is not an easy task in general to show that dynamics generated by 
physically-relevant equations preserve the convex set 

Ad := {q e Symo((i) : < XiiQ) < 1 - ^ for i = 1, d] , 

i.e. that Q{x, t) G A^^ for almost every x and t > 0. The following approach of [3], which 
has its roots in the paper of Katriel et al. [13], allows one to treat the modelling 
issue in both statics and dynamics in the same manner. 

1.4. The Ball-Majumdar Singular Potential. Recently, Ball and Majum- 
dar [3] proposed a qualitatively-similar continuum theory, for which the authors effec- 
tively 'build in' physicality of the Q-tensor. 

We now briefly outline their construction of a singular map -0 : SymQ((i) — )• MU{oo} 
and then introduce the problem in statics which motivates our present problem in 
dynamics. Although the construction of ip in [3] is performed in dimension d = 3 alone, 
it may be generalised in a straightforward manner to the case d = 2. 

1.5. The Roots of ip in Maier-Saupe Theory. As a starting point, consider 
the spatially-homogeneous Maier-Saupe mean field theory, 

ImsIp] '■=^ [ Pi^) logPl'^) duj+'^ [ [ (\- {u ■ a)A p{uj)p{a) duda, 
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for suitable probability density functions p : S'^^^ — > [0, 1] satisfying the symmetry 
property p{uj) = p{—uj). Moreover, 6 represents temperature and the constant k encodes 
to some extent information on molecular interactions. A routine calculation reveals that 
/ms has the form 

^21 



where Q is the Q-tensor corresponding to the probability density p. 

For a given physical Q-tensor Q G Symo((i), the authors consider the following 
natural entropy minimisation problem associated with /mSj namely 

min / p{uj) log p{uj) duj , (3) 

where 

Aq := \p -.8'^''^ [0,1] : [ p{u})doj = l and / (uj0U}-\l)p{uj)duj = Q 

In [3], it is shown that this problem possesses a unique minimising density p* in the 
class Aq, given explicitly by 

, exp(^ti^i^') 

Z(/ii, ...,Pd) 

where pi, ...jPd are Lagrange multipliers associated with the constraint that Q be the 
Q-tensor of the density, and Z{pi, ...,pd) is a normalisation factor which ensures p* is 
of unit mass. By uniqueness of minimisers of (3), one can then construct a related 
singular map : Symo((i) ^ M U {oo}. 



f 11 

/ p{uj)\og p{uj) duj if - - < \j{Q) <l - 



mm 



oo otherwise. 



Such a map establishes a framework in which non-physical Q-tensors are essentially 
forbidden. We henceforth denote the effective domain of V' by D(^) := {A £ Symo(d) : 
ip{A) < oo}. The following proposition records important properties of the potential ip. 

Proposition 1.1. The map ip : SymQ((i) — >■ MU{oo} has the following properties: 



Analytic Properties 



(PI) The map ijj is smooth on its effective domain, i.e. V ^ C°°{D{'ip)); 
(P2) It is bounded from below, i.e. there exists t/'o > such that — V'o < "^i^) for all 
X G Symo(rf); 

(P3) It exhibits logarithmic blow-up as X ^ (?D('0) from the interior. 



Geometric Property 



(P4) The map ip is convex on D{7p), which is itself a convex subset o/Symo(d). 



Algebraic Property 



(P5) The map ip is an isotropic function of d x d matrices, i.e. ip{RXR"'") = ip{X) 
for all R € SO(d), whenever X G Symo(d) is fixed. 
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Proof. We refer the reader to Ball and Majumdar [2] for details. 



□ 



Having constructed the map ■0, the authors replace the study of the Landau-de 
Gennes energy functional -ELdG with the following functional 



by replacing the Landau-de Gennes bulk energy density a(0)tr[Q^]/2 — 6tr[(5^]/3 + 
ctr[(5^]/4 with the singular map OipiQ) — Ktv[Q'^]/2 associated with the Maier-Saupe 
energy. The construction of ip in [3] is especially pleasing as whilst it resolves the above 
modelling issue for the static theory of nematics, the map ip is by no means confined 
to the study of statics. Using techniques from convex analysis (see Rockafellar [24] 
or Ekeland and Temam [10]) we study iIj in the setting of dynamics. 

It proves important to make the following dichotomy between physical and strictly 
physical tensor fields Q : — )• Symo((i). 

Definition 1.1 (Physicality). The field Q : Q, ^ Symo((i) is said to be physical 
on Vt if and only if —1/d < Xi{Q{x)) < 1 — 1/d for almost every x € 17 and i = 1, d. 

We compare this with the stronger notion of strict physicality. 

Definition 1.2 (Strict Physicality). The field Q : ^ Symo((i) is said to 
be strictly physical on $7 if and only if there exists 6 > sufficiently small such that 
5 — 1/d < Xi{Q{x)) < 1 — 1/d — 5 ior almost every x G and i = 1, d. 

The reader can check that if the map Q : Q ^ SyniQ^d) satisfies ip{Q) S for any 
given map satisfying (PI) to (P5) above, then it is physical. On the other hand, it is 
the case that Q is strictly physical if and only if 'ip{Q) € L°°. An issue similar to the 
eigenvalue constraint in liquid crystal theory also arises in the theory of elasticity (see 
[1]), where one is interested in demonstrating that the determinant of the deformation 
gradient tensor is bounded uniformly away from zero, i.e. det(Vy(x)) > 7 for some 
7 > 0. However, this problem is more challenging as constraints are being placed upon 
derivatives as opposed to the undifferentiated field variable. 

We now introduce the problem of study in [3] which motivates our present problem 
in dynamics from a modelling point of view. 

1.6. Motivation for the Problem in Dynamics. Consider the minimisation 
problem 



where the appropriate candidate maps in B have strictly physical trace on the boundary 
di}. One can show using a maximum principle approach that minimisers are strictly 
physical throughout the domain Q. We wish to answer a question which is similar in 
spirit for the case of dynamics. Supposing that the Q-tensor field evolves under system 
(S) above, we ask the following: 



Ebm[Q] :=^(^|Vg|2 + 0V'(Q)-^tr [Q']) dx 
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If initial data {Qo,uo) are of finite energy, namely 
£{Qo,uo) := J^^ (^|noP + ^|VQo|^ + ^^(Qo) - |tr [qI]^ dx < oo, 

is it the case that weak solutions Q{-,t) are strictly physical for t > 0, i.e. 
i^{Qi;t))€L^? 

This is one of the central questions of this paper and we answer it in the affirmative. 
Knowledge that ip{Q{-,t)) € L°° allows us to infer Q{-,t) E L°° by boundedness of the 
related eigenvalue maps. The 'strict physicality' property of solutions is an attractive 
feature from the point of view of regularity theory. Using the fact that the range of 
weak solutions Q then belongs to a compact subset of Symo((i), we are able to prove 
higher regularity of solutions that persists for all time t > in dimension 2. 

The unmodified version of system (S), taken from Beris and Edwards [4], has 
also been studied in the physics community by, among others, Yeomans et al. [8,26] 
and also in the mathematical community by Paicu and Zarnescu [20,21]. 

Remark 1.1. One might consider the above boxed question to be natural, as it 
is characteristic of evolution equations of parabolic type that solutions with initial 
data are instantaneously in L°°: for more details on such topics, one could consult 
Ladyzhenskaya, Solonnikov and Ural'ceva [14]. Although one may not be able 
to use the evolution equation for the tensor field to close an equation for the quantity 
ip{Q), one can nevertheless derive a parabolic inequation satisfied by ip{Q), from which 
a suitable comparison principle argument yields strict physicality of weak solutions Q. 



2. The Class of Models 

In what follows, the spatial dimension d will be either 2 or 3 and we employ the 
Einstein summation convention over repeated indices i,j,k,£,m and n with range in 
{l,...,(i} throughout. Moreover, (A) := A — d~^tT:[A]I denotes the trace-free part of 
any matrix A G M'^^'^. For strictly positive constants A,T, L,d, k and u, we search for 
maps Q : I'^ X (0,oo) ^ Symo(d) and u : I'^ x (0,oo) ^ R*^ on the spatial domain 
I'^ := [— A7r/2, A7r/2]^ which satisfy the coupled system 

^(x, t) + {u{x, t) ■ V) Q(x, t) - S{Q{x, t), Vn(x, t)) 



(S) 



= r (^LAQ{x,t) - (^|i(g(x,t))) +KQ{x,t: 
— (x, t) + {u{x, t) ■ V) u{x, t) = i^Au{x, t) - Vp{x, t) + div (r(x, t) + <t(x, t)) 



V • u{x,t) = 0, 

in the distributional sense. Solutions are subject to the periodic boundary conditions on 
the domain /"^ and evolve from given initial data Qq : I'^ ^ SymQ(d) and uq : I'^ ^ 
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to which they converge in an appropriate topology (strong-in-norm or weakly) as t tends 
to from above. 

The term S expresses to what extent the flow u locally 'twists' and 'stretches' the 
order parameter Q, and is given by 

S{Q, Vn) := {Do + ^D) (q + i/) - (q + i/) {Do - ^D) - 2^ (q + ^/) tr [QVu] , 

(4) 

where ^ € M is a rotational parameter whose value will be of some significance in the 
course of our analysis. The tensors Do and D defined by 

Do := \ ((Vn) - {^uf) and D := ^ ((Vtx) + {Vuf) 

are the anti-symmetric and symmetric parts of the velocity gradient tensor Vu, respec- 
tively. The stress tensors r and a are given component-wise by 



dQ dQ 

dxi dxj 



(5) 



and 

c^ij := QikHkj - HikQkj, (6) 
where H is defined for notational simplicity to be 

H :=LAQ-e(^{Q)) +kQ. 

Moreover, the singular map ip belongs to the non-empty class of all maps on SymQ((i) 
satisfying properties (PI) to (P5) above. 

Remark 2.1. The right-hand side of the Q-equation is the formal gradient of the 
energy functional Ebm which replaces the usual contribution from the more commonly- 
adopted Landau-de Gennes energy ii^LdG- On thermodynamic grounds, we justify such a 
replacement since both bulk potentials are smooth on their respective effective domains, 
and are qualitatively similar to one another, in the sense that both can describe a first- 
order nematic-isotropic phase transition (see Section 4 of [3] for more on this point) 
and possess the same material symmetry, viz. (P5). 

2.1. Statement of Main Results. The following two theorems contain the main 
results of this paper. 

Theorem 2.1. For initial data {Qo,uo) G x L^^^ of finite energy 
£{Qo,uo) = J^^ (^|uoP + ^iVQoP + dij{Qo) - |tr [qI]) dx < oo, 

there exist maps Q G L'^^{{0, oo); H^) n Lf„^((0, oo); F^) and u € L'^^{{<d,oo);L\^) n 
-^k)c((^' '^)' -^div) ^/^^c/i satisfy the coupled system (S) in the distributional sense. The 
map Q is also strictly physical for positive time, i.e. 

ip{Q{-,t)) G L°° for almost every t > 0. 
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Thus, if one assumes £{Qq, uq) is finite, one can infer the strict physicahty of weak 
solutions Q{-,t) for t > 0. If one endows the initial data with higher regularity and 
insists on Qq being strictly physical, the following holds true. 

Theorem 2.2 (dimension d = 2). If uq e //^iv ^^'^ Qo G ^'^ with ijjiQo) G 
then there exist maps Q £ L'^{0,T; H"^) D L'^{0,T; H^) and u G L°°{0,T; H^^) n 
L^(0, T; iJjjy) for any T > which satisfy (S) in the sense of distributions. 

2.2. Regularisation of ip. For the purposes of building a weak solution to the 
system (S), it is inconvenient to work directly with the singular map 'i/'- We work instead 
with a regularised map ipN parameterised by the mollification index A^ = l,2,3,..., such 
that ip can be recovered in the limit N ^ oo. 

Firstly, for J = 1,2,3,... we define ipj : SymQ((i) — )• M to be the Yosida-Moreau 
regularisation of namely 

MQ) ■■= , min (J\A-Q\' + ^PiA)) . (7) 

Secondly, for fixed J and for any K = 1,2,3,... we define ipj^K to be the standard 
mollification of the map ipj, namely 

^J,KiQ) := K''' [ MKiQ - R))HR) dR, (8) 

where ^ G C^(M'^^'^, M+) has the unit mass property J^d^^d ^{R)dR = 1. Finally, we 
define tpN '■= fpN,N for each > 1. We now quote without proof a number of properties 
oi ipN, taken from Feireisl et Al. [11], which are of use to us in the construction of 
weak solutions. 

Theorem 2.3. For each N > 1, the regularisation tp^ of the Ball-Majumdar 
potential has the following properties: 

(Ml) The map ipN is both C°° and convex on W^^^; 

(M2) It is bounded from below, i.e. —-00 < iPn{X) for all X G M^^"' and for all 

N >1, where "00 > is the same constant appearing in (P2); 
(M3) Vtv < i^N+i on M.'^'"^ for N > 1; 

(M4) -i/^AT —^tpin LJ^^(D('i/;)) as N ^ oo, and ip^- is uniformly divergent on SymQ{d)\ 
□ (-0) as N ^ oo; 

(M5) ^ ^ |i zn LZiDm as iV ^ oo; 

(M6) The regularised map ipN satisfies 



dQ 



< clf\x\ + cff 



for all X G M''^'^ and positive constants c^y and C\q which depend on the molli- 
fication parameter N > 1. 
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3. A Priori Estimates for the System (S) 

As is customary, we precede in a formal manner to obtain useful a priori estimates 
associated with the system (S). The following manipulations hold only for smooth maps 
Q and u, however the resulting energy estimates associated with (S) do indeed hold for 
suitable approximants Q^^^ and u^^^ in Section 4 to come. Moreover, for convenience 
we replace ip with a convex map '0* : Symo((i) — [—^'0,00) satisfying properties (Ml) 
to (M6) above. Consequently, we shall have no uniform information on maps Q^^^ 
in space. We recover such information in the limit: see Section 4.3 below. 

3.1. An Energy Identity for the System (S). Firstly, considering the evolution 
equation for the Q-tensor field, 

^ + (n • V)Q - SiQ, Vn) = T (lAQ - 9 {^{Q)) + 

and taking inner products throughout against —H, where 

H = LAQ-e{^{Q))+KQ, 

one finds after integrating by parts that 

K^IIVQII^ + MQ) - f IIQIlO +r LAQ-e {^{Q)) + ^q"' 



dQ 



2 



L J^^{u-V)Q : AQdx + J^^iiVoQ-QoJo) ■ {J^ (Q)) ' O) 



-L J ^ {uJoQ - Qojo) : AQdx- J ^ [oji [q + i/) + (q + ^l) oji) ■ AQ dx 
+ 2Li J ^ tr [QVn] (q + i/) : AQ dx 



Secondly, considering the evolution of the velocity profile 

du 

— + {u ■ V)u — vAu + Vp = div (r + o") 

and this time taking inner products in L^j^ throughout the equation against one 
obtains in a similar way by parts 

^|(lHli) + H|v.||i 
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L j diQjidkQejdkUi dx — L J (QAQ — AQ Q) : Vudx 



2Li j^^ tr [QAQ] [q - i/) : Vu dx 



7-5- 



+ 2e/^tr[Q(^( 



(Q) 



/ : Vn (ix 



7;- 



{Q)) -kQ]Q] -.Vudx. 



(10) 



J 



Let us write £ : x L^j^ 



^ M U {00} to denote the functional 
£{Q,u) ■.= ^\\VQ\\l + e J^^^l;,{Q)dx-^\\Q\\l + ^\\u\\l 
By adding the contributions of (9) and (10) together and noting the cancellations 



Tj +Tj = for j = 1, 6 together with the null terms 2 = ^ = 0, one discovers 



d£ 
'dt 



-r 



LAQ-e(^^{Q)) +kQ 



i^ll Vn||2, 



(11) 



which implies that the functional £ is non-increasing along solution trajectories {Q{t), u{t)) 
in X Ll,^. 

At this point, we have paid a price for the cancellations of higher derivative terms 
(like Ti' and , for example) which would otherwise complicate our subsequent anal- 
ysis: the functional £ is not of one sign due to the negative contribution of the 2-norm 
of the tensor field. This would cause a problem in the sequel if we wish to use £ to 
obtain bounds on Q and u and their distributional gradients in natural function spaces 
which are uniform in approximation parameters. 

This issue, however, is remedied as follows. Suppose now that the smooth maps 
Q = and u = u^"> depend on an approximation parameter /3 G {1,2,3,...}, but 
that the initial data Qo S and uq € are independent of /?. Equality (11) above 
implies the simple inequality 

|(f(g('^)(t),^.(/^)(t)))<o, 

from which we deduce by integration in time, 



12 

L 



< ol|VQo||2 + ^ / MQo{x))dx-^\\Q4l + \\uo\\l 

I J Id I 

Now, noting that satisfies the property (M2), which impUes that 



(12) 



Jid 

and if we have in addition that 

sup||Q(«(-,t)i|2 <co, 



(13) 



where cq > is some constant independent of /3, we may deduce from (12) using such 
bounds that 



sup 

t 



L 



l|VQ 



i||n(^)(,t)||i)<Co, 



(14) 



where 



L 



VQolli + \\\uo\\l + f (^0 + MQo{x)) dx + co> 0. 

Z J Jd v ^ ' 



>0 



Thus, the existence of the functional £ along with the additional stipulation that 
{Q^'^^lflLi be uniformly bounded in Lf^L^. allow us to infer that 



{Q^^^lfci is uniformly bounded in Lf'H^ 



OO Tjl 



(15) 



and 



is uniformly bounded in LfL^^. (16) 

We can use these uniform bounds to glean yet more information from the identity (11). 
Performing an expansion of the right-hand side, one finds 



d 



(^(Q(^)(t),n(^)(t))) +rL2||AQ(^)||2 + r02 (^(Q(/^))) 



+ Tn\\Q 



|Vt.(^)||2 + 2rL/^||VQ(^)||2 + 2rL0 / AQ(« : (^{Q'^^^)) dx 

J Id \ oQ I 



+ 2TeK 



Id 



\dQ 



{Q^Py)) -.Q^^Ux. 



(17) 



Now, by integration by parts, we obtain 



AQ(^):^(g(^))dx 



d 

E 

m=l 



Since the map ip^ is convex, its Hessian is positive definite: see, for instance, ROCK- 
AFELLAR ([24] Chapter 4, Theorem 4.5). Using the convexity of ■0) one may deduce 
from the above that 

2rL6l 



AQW:^iQ(f^^)dx<0. 



By Young's inequality, one may also obtain the simple estimate 
lidXdQ 



2reK 



Q^'^ dx<^ 



/^^*(Q(^); 



\dQ 



+ 2Tk^\\Q 



Using these observations, integrating in time across equality (17) yields: 



{(3^'^^}^=i is uniformly bounded in L^H^ and LfH'^, 
{u^^^'^^i is uniformly bounded in L'^L'^i^ and L^H^^^, 



(18) 
(19) 
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and 



\dQ 



is uniformly bounded in Lfll. (20) 



13=1 

Thus, the moral of this section is that if we can establish identity (11) rigorously, along 
with L^L^-bounds on {Q^^^}f^i and the 'convexity' inequality 

2rLe I AQ(''):^(Q(''))dx<0, (21) 

we are immediately in the familiar setting in which we aim to employ weak compactness 
arguments to identify candidate solutions for a limiting system. 

3.2. An Estimate for Higher Regularity of Solutions. The property of strict 
physicality becomes of importance when investigating higher regularity of weak solu- 
tions of (S). In particular, strict physicality of weak solutions Q{-,t) for t > implies 
that ip{Q{-,t)) is as distributionally differentiable as Q{-,t), since the range of the tensor 
field belongs to a fixed compact subset of Symo(d). 

In what follows, we only consider the case of two spatial dimensions. By noting 
convenient cancellations in the system (S), knowledge that ip{Q{-,t)) E for t > al- 
lows us to prove Q G L'^{0,T; H^)r\L^{0,T; H^) and u G L'^{0,T; Hl^)r\L^{0,T; Hj-J 
for any T > whenever initial data {Qo,uq) lie in the smaller class x H^-^^ with Qq 
satsifying ip{Qo) G L°°. 

Our remarks here are once more merely formal, but shall be made rigorous in 
Section 5. Consider the auxiliary functional T : x H^-^^ — )• [0, oo) defined by 

F{Q,u):=-( \Vu\^dx + -f lAQpdx. 

2 Jjd 2 J Id 

We assume once again that the maps and u^'') are smooth, with the additional 
stipulation that Q^^\-,t) be strictly physical on for t > 0. By considering the time 
derivative of T'^'^Ht) := T{Q^^\-,t),u'^^\-,t)) and noting the identity 

L f Auf^ ( dx-L f AQ(/5) : A ((u(« • V)Q^m dx 
Jjd dxj \ oxi oxj I J Id ^ ' 



one can show that 

dT^ 
dt 



+ v [ \An^('^\Ux + TL^ f |VAQ(^)|2dx 

Jid Jid 

= [ {u'-^^ ■V)u^'^^ ■ Au^'^Ux + TLk [ iVQ^^^'^dx 

Jid Jid 

+2l[ ^(AQ(^))..(uW.V)^dx-rL9 [ AQ^ : a((^{QW))) dx. 
Jiddxk oxk Jid \\oQ I J 

Now, under the assumption that Q'^^^ is strictly physical, it can be shown by means of 
Ladyzhenskaya's inequality that 

^{t)<CoT^^\t)^ + Ci, (22) 
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where Co,Ci > are independent of the approximation parameter /3. Since IS 
uniformly bounded in L^(0, T) by (18) and (19), an apphcation of Gronwall's inequaUty 
yields that the approximants Q^^^ and u^^^ are uniformly bounded in L°°{0,T; H^) n 
L2(0, T; H^) and L°°(0, T; iJ^j^) n ^^(0, T; Hj.^) for any T > 0. Such improved uniform 
bounds give us better compactness and convergence properties when passing to weak 
solutions in the limit /3 — )• oo. 

Now that we have made our initial remarks regarding the formal structure of the 
system (S), we embark upon a construction of weak solutions. 

4. Existence of Weak Solutions 

We posit a family of coupled systems which may be considered as approximants to 
the coupled system (S). Weak compactness arguments identify candidate solutions to 
the 'limit' system (S), which satisfy the coupled system distributionally in the limit by 
strong compactness of the set of approximate solutions. 

Stage I to Stage IV below should be considered as a schematic for the proof of 
Theorem 2.1 above. The rest of the details may be sourced from texts such as Lions 
[17] or Robinson [23]. 

Remark 4.1. For simplicity, we only consider the case ^ = as the approximation 
scheme we adopt in this case is compatible with the strict physicality argument in Stage 
III below. We pass comment on the reason our comparison principle argument is not 
suitable for the non co-rotational case in due course. 

Stage I: Regularised Potential and Finite-dimensional Velocity Field 

Let us set out the basic objects which which we employ in this stage. For G M, 
consider the Stokes operator eigenfunction problem on /"^ = [— A7r/2, A7r/2]'^ given by 

— z^A(/)(x) + Vp(x) = iJ-(p{x), 
V • (t){x) = 0, 

supplemented with periodic boundary conditions on I'^. It is known (see, for ex- 
ample Temam [25]) that there exists a countably-infinite family of eigenfunctions 
(pij with corresponding eigenvalues fiij := 47r^z^jj|^A^ of multiplicity m{i) for in- 
dices i,j > 1. We enumerate this countable family of eigenfunctions and eigenval- 
ues by {(j)m}m=i ^'^d {//m}m=i) respectively. It is also well known that the fam- 
ily {4>rn}m=i Constitutes an orthonormal basis for L^^^ and an orthogonal basis for 
i/j;^. Define Hm ■= span{(/)j}|£^ and the associated orthogonal projection operator 

■■ ^div ^ by ¥mu := J^kU i^, (pk)2 0fc for u £ L^j^. 

We now state the main result of this stage. 

Proposition 4.2. For any uq G L^;^, Qq G satisfying ipiQo) G and any 
fixed positive integer M > 1, there exist maps Q^^'^^ G C°^{I'^ x (0, oo); SymQ(d)) and 
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y^{M,N) g C'-'^((0, oo); //m) satisfying the system (Sm) given by 
V • n(A^'^) = 

pointwise in Symo((i) anrf M*^ respectively, for all x £ I'^ and t G (0, oo). 

Proof. The result follows from a rather involved Banach and Schauder fixed-point 
argument, whose proof we omit. The basic form of the argument can be found in Lin 
AND Liu [16]. □ 



Stage II: Passing to the Limit — > oo 

We now aim to show that our approximants Q^^^'^^ and n^*^'^-* are uniformly 
bounded in in natural function spaces associated with the a priori estimates of Sec- 
tion 3. By compactness methods (such as those of Banach- Alaoglu or Aubin-Lions), we 
may then identify candidate maps which are solutions of a limiting system of equations 
(Sm) as the parameter N oo. Throughout this stage, M > 1 is chosen arbitrarily at 
the start and remains fixed. 

4.1. Uniform Bounds. In order to carry out the program sketched in Section 3, 
we need only demonstrate that is uniformly bounded in Lf^L^, since our 

approximants are sufficiently regular that both the identity (11) holds and 

d 



jd 



E 

m=l 



^{M,N) 



by an application of classical integration by parts, since iI^m '■ Symg((i) — )• [— 1^0,00) is 
smooth and convex by property (Ml). 

We consider the evolution equation for and recast it in distributional form, 

namely 



Id \ dt 



+ (n(^'^) • V)Q(*-^'^) - 5(g(^^'^), Vtx(*^'^)) ]:xdx 



r (laq(*^'^) - 9 



ijd 

for any x £ ^t^x- Choosing x to be the element 



dQ 



+ kQ 



{M,N) 



: xdx, 



(23) 



one may derive the equality 
d / K 



dt\2 
TL9 I 



Jd OQ 



(M,N) ||2 
Il2- 



(24) 
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Now, since property (M2) of the mollified map -i/^Tv implies that 

one may deduce from (24) that 

d_ 
dt 

and since by property (M3) 



^Af(Qo) dx <e ^(Qo) dx, 

J Id 

we find that is uniformly bounded in Lf^Ll.. Applying the reasoning of 

Section 3, we quickly deduce that 

|g(M,iV)|oo^^ is uniformly bounded in L'^hI and lIhI,, (25) 
|^(M,Af)|oo^^ is uniformly bounded in L'^L\^ and LIhI^^ (26) 

and 

'^(Q(*^'^)))|°° is uniformly bounded in L^L^. (27) 

/ J N=i 

By Banach-Alaoglu compactness, we may extract subsequences of (25), (26) and (27) 
which converge weakly to limit points Q(*^), n^^^) and Y^^^^ in LjHl, 1^1%^ and L^Ll, 
respectively. 

We may in fact extract more information from the boundedness property of {n^*^'^) }?^=i 
above. Firstly, by orthonormality of the family {(t)k}'kLi i'^ -^div' 2-norm of u^^'^\-,t) 
is realised as 

||^(^W(.,i)||2=;^|4M,iV)(,)|2^ 
r7i=l 

a quantity which we know to be uniformly bounded in N . Using the fact that one may 
construct a norm which is equivalent to the standard norm on -ff^j^ using fractional 
powers of the Stokes operator, one may use the above observation, together with (26), 
to show that 

|^(M,Af)|oo^^ is uniformly bounded in C7([0, T]; //J;^) for any s G R. (28) 

The uniform bounds of (25) and (28) and allow us to show in turn directly from the 
Qi^^^)-eqnQX.iou that 

i ^® uniformly bounded in LIL^ 

and due to the higher spatial derivatives of Q^^^'^^ in the forcing term of the velocity 
field equation, the weaker statement that 

oo 

is uniformly bounded in L1H7-^ 



I dt 



t div ■ 

N=l 



One may verify that the weak limits of {^tQ^^-^'^^j^^i in LIlI and in 
LjH^^ are dtQ^'^^ and dM^\ respectively. 
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4.2. Strengthening Convergence. For the purposes of passing to a suitable 
limiting equation, we must now strengthen our notions of convergence of approximants. 
For the case of the tensor field, by virtue of the chain of compact embeddings CC 

CC L^, Aubin-Lions compactness allows us to strengthen the weak convergence of 
Qi^'i^^) to its limit Q'^'^^^ to strong convergence in LjHl. For the case of the velocity 
field, since the time derivatives of the co-ordinates of u^^^'^^ satisfy 

rt dc^^^^) 

c(^^^)(t) = {uo,cPk)2ek + ^^is)ds for < t < T, 

one may show that {c^^'^'^^ j^^^^ constitutes an equi-continuous family of maps in the 
space C([0, r]; R*^). By Arzela-Ascoli compactness, we deduce that c^*^'^-* converges 
strongly to some c*^^^) in C([0, T]; M^^), following which one may deduce that n*^^^) is 
equal to c^^^Vfc- Furthermore, converges strongly to u^^^ in C{[0,T]; H^-^^) for 

any s G M. In particular, we have that u^^) in C{[0,T] x I'^). 

4.3. The Range of the Map {x,t) i— > Q^^''^\x,t). In order recover pointwise 
information on the weak limit G L'^L'^, we must show that Q^^^\x,t) € 
almost everywhere. To this end, it is sufficient to show that 

ip{Q'-^'^\x,t))dx < oo 

Jd 

for almost every < t < T, as 'ip{Q^^^\-,t)) being integrable implies that Q^'^-^\x,t) G 
D{ip) almost everywhere. 

We fix No > 1, and suppose N > Nq. Since we know that ^ Q^^^ in L^i?^ 

by section 4.2, it follows by the smoothness property (Ml) of the mollified map Va^o 
that 'ipNoiQ^^^'^\xyt)) ~^ '4'No{Q^^\x,t)) almost everywhere as ^ cxd. Recalling 
the uniform bound 

" i^N{Q^^'''Hx,t))dx<Co, 



lid 

for some constant Cq > independent of both N and time, by Fatou's lemma and the 
monotonicity property (M3) we may deduce that 

^No{Q'^^'Hx,t))dx<Co, (29) 

Ijd 

almost everywhere in time. An application of the monotone convergence theorem finally 
allows us to deduce that Q^^^\x,t) G D('0) almost everywhere. 

4.4. The Limiting System as N ^ oo. We now consider passing to a limiting 
equation which the weak limits Q^^'^^ and u^^^^ satisfy. It follows that since Q^^'^\x, t) 
tends to an element of D^ip) almost everywhere, 

^(Q(^'^)(x,t)) ^ ^(Q(^^)(x,t)) almost everywhere on /'^x(0,r), 

by property (M5). By a dominated convergence argument, one may conclude that 
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Finally, passing to the limit N ^ (yo we find that the maps 

q(m) g L^o^L^ n 

LfHl n LfH^ and u^^^ G -Lf Ljiv ^?-f^div satisfy the coupled system 



(Sm)< 



+ (u(^) • V)n(*^) - z.Au(*^) + Vp = div (r(Q(^)) + 



dt 

V • tx(^) = 



in the distributional sense. One may also swiftly verify by standard methods that 
G C([0,r];L2) and G C7([0, T]; (L^^J,). 

Stage III: Strict Physicality of Approximants 

In the previous stage, we were able to demonstrate that Q^^\x,t) G D('0) almost 
everywhere, which was a necessary step before passing to the limit system (Sm) above. 
We can in fact show more, namely that there exists 5 > which is independent of 
M > 1 such that 

<5 - ^ < Xj{Q''^\x,t)) < 1-^-6, for almost every and j = 1, ...,d, (30) 

i.e. approximate solutions Q^^^ are strictly physical for almost every t > 0. We shall 
obtain strict physicality of approximants Q^^^ (and, in turn, of weak solutions in the 
limit as M — 7- oo) by a suitable comparison principle argument. We apply the maximum 
principle at this stage in the construction, as weak solutions lack sufficient regularity 
for these methods to be applied. 

4.5. A Parabolic Inequality for ipNiQ^^'^^)- We begin at the level of regularity 
obtained at the end of Stage I. We once more recast the evolution equation for 
in distributional form, namely 

£ J^^ (^^^ + (^^"^'"^^ • V)Q(*^'^) - 5(Q(^^'^), Vn(^'^))) : xdxdt 



= r £ J^^ (lAQW -0(^(Q(^'^))) +AcQ(^'^)) : xdxdt, 
for any x G -^t-^^- Making the following choice of test function 

X=(^(Q(^'^)))¥' for any <^ G C7-(/'^ x [0, T], 
one may quickly deduce that the inequality 

I (V'iv(Q(^'^))) + (n(^'^) • V)g(^'^) - rLAViv(Q(^'^)) < ^tr [(Q(^'^))'] (31) 

holds pointwise in M for all {x,t) G I*^ x (0, T). Using a suitable comparison function, 
we now demonstrate that it is possible to control i^NiQ^^'^'') in uniformly in M 
once we have passed to the limit N ^ oo. This uniform control ultimately allows us to 
carry information on ip(Q^^^) through to the weak limit il^{Q). 
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Consider now the maps G^^'^) , i7 : I'^ x (0,T) R, where G^*^'^) solves the 
homogeneous problem 

QQiPI,N) 



+ • V)G(^^'^) - rLAG(^'^) = 0, 



N) 



with mean zero initial data, and //(^^'^) solves the inhomogeneous problem 
(P^) 



p,u{M,N) Yk"^ r o 

+ (u(M,N) . y^H{M,N) _ riAi7(^^'^) = — tr (q(^^'^))' 



with constant initial data. Both problems (P]v) and (P^) are supplemented with 
periodic boundary conditions on Defining the map 

K(M,iV) to be the difference 

V'wlQ^*^'^^) - G(^'^) - H^^^^\ one quickly sees that satisfies the parabohc 

inequality 

^^^^ + (n(^^'^) • V)E:(^^'^) - rLAi^(^^'^) < 0, (32) 

pointwise on I'^ X (0,T) with i^(*^'^)(-,0) = 0. By the classical parabolic maximum 
principle (see, for example, Pucci and Serrin [22]), we find that (32) above implies 
that 

K(^'^)(x,t) < on /^x(0,r) 
^ Vjv(Q('^''^H2;,i)) <G'(*'''^)(x,t) + F(^^'^)(x,t) on /'^x(0,r). 

Thus, if we can obtain bounds independent of M on G*-^^'^-* and //(*^'^) in in the 
limit as ^ oo, we can infer strict physicality of both approximants Q^*^^ and, in 
turn, weak solutions Q as M ^ oo. We now perform such an analysis on these two 
comparison functions. 

4.6. Analysis of the Comparison Function G^^'^\ By a standard construc- 
tion, one can show that problem (Pjy) has a unique solution which is classically smooth 
for t > 0. Although we know that the approximants ipNiQo) lie in by property (M6), 
we only know V'(Qo) to be in L^. Thus, to gain uniform control on G^*^'^^ in L°° in 
the limit as — >■ oo, we require the following — >■ L°° estimate from Constantin 

ET Al. [6]. 

Lemma 4.3. Let v be a smooth, spatially-periodic divergence-free velocity field v, 
and let 7 > 0. Suppose that g evolves under the associated advection- diffusion equation 
on the two- or three-dimensional torus, namely 

^ + {v.V)g- TLAg = 0, 
g{-,0) =50 G L\ 
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where go is of zero mean over I'^. There exists a constant C = C(7) > which is 
independent of v such that g satisfies 



Il5'(-;' 



< — — l|5o||i 



t2 



+7 



(33) 



for t > 0. 

Applying the estimate (33) to the solution of (P]v)) one can show that 



< 



C 



t2 



1+7 



iPn{Qo 



1 



(Avr)'^ Jjd 



^P{Qoiy))dy 



for all x G I*^ and t > 0, which together with property (M3) of the regularised potential 
yields 



< 



C 



t2 



+7 



IIV'(Qo)l|i, 



(34) 



where the resulting bound is clearly independent of both M and A^. 



4.7. Analysis of the Comparison Function H^^^'^\ We now compare the 
smooth solutions of problem (P^) with those of the 'limiting' problem 



(P^ 



dt 



(Avr) 



Id 



26 

V'(Qo) dx 



{Q 



(M)^ 



Since the initial datum is constant, Q^*^) G L'^H^nLfH^ with Q^^^\-,t) G L°°, and the 
vector field u^^^^ is also smooth in space and time, it follows that F^^^) G L'^H^nL'^H^. 
In particular, we have by the Sobolev embedding theorem that H^^'^\-,t) G Cpgj.(/'^) 
for < t < r. If we denote the difference ifC^'^) - H^-^'^ by H^^^^\ a can be checked 
that it satisfies the equality 



dt 



{M,N) 



FLAH 



{M,N) 



-tr 



29 

l^iM^N) _ n(^^)) . V) iJ^^^'^) 
f{M,N) 



{Q 



{M,N) 



(M) 



in for each t > 0. Multiplying throughout by H^'"'" ' and integrating over the spatial 
domain / , one may use the fact that is uniformly bounded in LfLl to 

derive the inequality 



i^|i;^(M,iV) 



\l + FL\\WH^'''''\;t)\\l<C\\Q^^''''^ 



{■,t)-Q^''\;t)\\ 



+C\\H 



{M,N) , 



.t)||i + C (max||n(*^'^)(,.) - uW(,.)|U ) ||V//(^^ 



which holds for < t < T. Finally, an application of Gronwall's inequality along with 
the results of section 4.2 yield the convergence result 

hm ||//(^^'^)(-,t)-if(*^)(-,t)||2 =0, 

from which we deduce that 

H^^^'^\x,t) H^^^'^\x,t) almost everywhere on /'^x(0,r) as A^ ^ oo. (35) 
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At this point, we can say the following. Since 'i/'MlQ^^'^H^; 0) ^ G^^^'^\x,t) + 
H^^''^'^\x, t), using property (M4) of the regularised potential and the fact that Q^^-^'^\x, t) 
Qi'^'0(^x,t) € almost everywhere as ^ oo, we may deduce that 

V'^(g(*^'^)(^,t)) < ^11^^(^0)111 + 

= lim MQ^''^^\x,t)) < -^\\^l;iQo)\\,+H^'''\x,t) 

for almost every (x, t) el'^x (0, T) by results (34) and (35) above. Therefore, if we can 
find a bound on //(*^) in L°° which is independent of M, we have demonstrated our 
earlier claim that tp{Q^^'^\-,t)) may be controlled in uniformly in M. 

4.8. Uniform Bounds in on H^^^\ As solutions of problem (P^) are suffi- 
ciently regular, multiplying throughout the equation by for p > 2 
and integrating over I'^, we find 



dx 



'29 Jjd 



whence 

1 d 

p dt 



p Jl'l 



29 



tr 



(Q 



1 d 



p dt 

|^(^')(-,t)||p<2F 



~29 



2\P 



(Avr) 



+ 1 



dx + (l 
1 



1 



e(i-i/p)T + 



29 



p 

FW(,t)||^ 



i7W(.,t)||^ 



{p — l)p 



Finally, using the fact that hes in LP for every p S [1, oo) U {oo}, by taking the 

limit p — > cx) in the above inequality, we deduce that solutions of (P^) satisfy 



'il^{Qo)dx 



29 



(36) 



and so the sequence {H^^^^}^^^-^ is indeed uniformly bounded in L'^L'^ as claimed. 
Piecing together the remarks of this stage, we conclude that 



HQ^''\x,t))<^mQo)h + 



jd 



ipiQo) dx 



29 



(37) 



whence Q^^'^\-,t) is strictly physical for t > 0. Such a property will automatically be 
inherited by weak solutions Q if we can show Q'^^'^\x,t) — > Q{x,t) almost everywhere 
as M — )• oo. 



Remark 4.1. Let us comment briefly on the non co-rotational case € M \ {0}. 
If one follows through the scheme of Stage III, one finds that the relevant comparison 
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function H in Section 4.7 should satisfy 



where co(0) = 0. Our construction of weak solutions does not provide us with VFJj^- 
bounds on approximants u^^'^ uniform in M which we need to carry out the program 
outlined above. It is for this reason we restrict our study of strict physicality of weak 
solutions to the case = 0. 

Stage IV: Passing to the Limit M — > oo 

At this point, compared with Stage II we have much less to do in order to identify 
candidate maps for weak solutions. Since we have shown in the previous stages that 
Q^^'^\x,t) G D('i/') almost everywhere, it follows that Q(*^) g Lf^L^ and so is auto- 
matically in L'^L'^. Furthermore, our maps g(^^) and 

^(Af) sufficiently regular that 
the identity (11) holds. It is also now straightforward to verify that 

AQ^''^:^{Q^''^)dx (38) 

is of positive sign. Since we know the image of the maps Q^^^^ belong to a compact 
subset of D(^) which is independent of M, we need not worry about distributional 
differentiability of dQip{Q^'^^). Once again, the reasoning of Section 3 allows us to infer 
that 

{Q^^^^}m=i is uniformly bounded in L^H^nLfH^, (39) 
is uniformly bounded in L^L%^ n L^H^^ (40) 



and also 



{|^(g(^^))| is uniformly bounded in ifll. 



(41) 



' M=l 

We shall denote the weak limit points of the sequences (40), (26) and (41) by Q, u and 
Y, respectively. 

Using the fact that the approximate tensor field equation holds in the strong sense 
in LfL"^, we may verify that 

{~^~} is uniformly bounded in LfHP^'^\ (42) 

where p{2) = and p{3) = —1. Similarly, considering the distributional form of the 
approximate velocity field equation, one may in turn check that 



is uniformly bounded in H^^. (43) 

M=l 

To strengthen the convergence of the approximants Q^^^) and u^^^^ to their weak limits, 
we need only apply once again the Aubin-Lions compactness lemma. 
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Finally, using the facts that ^ Q in ^H^, n^^) ^ n in LfL^-^ and that 

both Q^^'^\x,t) and Q{x,t) belong to a compact subset of D(^/^) almost everywhere, 
passing to the limit in system (Sa/) we can show that the system 

' dtQ + {u-V)Q-S{Q,Vu)=r(^LAQ-e(^^{Q)^ +kQ^ , 

^ dtu + (n • V)u + Vp = uAu + div (r + a) , 

V ■u = 

is satisfied distributionally hy Q e L^H^D L^Hl and u G L^L^-^ D if-f^div 

This completes the construction of weak solutions and, in turn, closes the proof of 
Theorem 2.1. 



5. Higher Regularity of Weak Solutions in Dimension 2 

As intimated in Section 3.2, if one places higher regularity conditions on the initial 
data for the system (S), strict physicality enables one to prove that weak solutions of 
Theorem 2.1 are in turn more regular. Furthermore, higher regularity of solutions 
allows one to infer that the limiting Q equation holds in the strong sense as an equality 
in ^Ll. 

We begin our enquiries at the end of Stage III, under the additional assumption 
that uq € H^[y and Qo G with tp{Qo) £ In particular, the reasoning of this 

stage (in particular, estimate (37)) implies that il){Q^^'^\-,t)) G L°° for t > 0. Moreover, 
by comparing approximants Q^^^^ with strong solutions i?^^'^^ of the problem 

^^^"^^ FLAijW = -(nW . V)Q(*^) +5(Q(*^),VnW) 



dt 



i?W(.,0) = Qo, 

one may deduce by uniqueness that Q(^^) eL^H^n LlHl for each M > 1. 
With this in mind, consider the 'higher-order' energy 

^W(t):=^(QW(.,t),n(*^)(.,t)) = l / |VnW(x,t)|2dx + - / | AQW(x, t)|2 dx. 

2 J [d 2 Jjd 

Using the identity 
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one can show that, at this level of regularity, the following energy identity holds: 

at Jjd Jjd 

= I • V)n(^^) • An(^^) dx +TLk ( IVQ^^'^^P dx 

Jid Jjd 
v ' 

/Ci: = 

+ 2l[ {^u^m.V)^QW,^^Qmdx-TLe I AQW:a((||;(QW))) dx. 
Jjd dxk oxk J Id \\oQ I J 

K.2: = 

By means of Ladyzhenskaya's inequality and the uniform bounds (39) and (40), one 
may derive the estimate 

X;i<^||AnW||i + C(||VnW||^ + l) 

and also 

X:2<£^||VAQW||2 + c(||VnW||^ + ||AQ(^)||9, 

for constants C > which are independent of M. Furthermore, since strict physicality 
implies the existence of a compact subset K C SymQ((i) such that Q^'^'-^\x,t) E K for 
all M > 1 and a.e. {x,t) G x [0,T], one has by property (PI) of tp that 



^ {Q(^^\x,t)) 



< Ck for k G {1,2,3}, 



for some constant independent of M. Such information leads one to deduce that 

^(t)<Co(^W(t))%Ci (44) 

for M-independent constants Ci > 0. Multiplying throughout (44) by the integration 
factor 

exp (^-Co J^^^^\s) ds 

acknowledging left-continuity of at and using the fact that J^(^) is uniformly 

bounded in (0, T) for any T > 0, an application of Gronwall's inequality allows one 
to deduce that 

{'3^^^^}m=i is uniformly bounded in L^H^nlfH^, (45) 

and 

{^^*^^}m=i is uniformly bounded in L^H^^nLfHl^. (46) 

Such uniform bounds give us more with which to work as we pass to limiting 
weak solutions of system (S). Aubin-Lions compactness guarantees us the existence of 
(relabeled) subsequences of {(5^^^''}m=i ^^'^ {u^^^^}m=i which are strongly convergent 
in L^H^ and L'fH^^^, respectively. In particular, one may show that {dtQ^^^}'^=i 
is uniformly bounded in L'^L'^, and further that the limiting Q equation holds in the 
strong sense as an equality in LIL^.. 
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